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Bose- Einstein coalescence of independent photons at the surface of a beam splitter is the physical 
process that allows linear optical quantum gates to be built. When distinct parametric down- 
conversion events are used as an independent photon source, distinguishability arises form the 
energy correlation of each photon with its twin. We derive upper bound for the entanglement which 
can be generated under these conditions. PACS: 03.67.Mn, 42.65.Lm, 42. 50. St. 

I. INTRODUCTION 

Interference of photons emitted by independent sources is one of the most intriguing features of the quantum 
theory of light and represents a remarkable departure from classical electromagnetism p], This phenomenon is 
easily understood if the radiation is treated in terms of creation and annihilation operators, rather than by invoking 
the superposition principle for light waves. This purely quantum effect represents the key feature for entanglement 
generation needed to build linear optical quantum gates [3J. These schemes rely on Hong-Ou-Mandel (HOM) effect, 
which is the observed coalescence of photons on the outputs of a beam splitter when they arrive simultaneously on 
the two input ports 0]. Since this effect is due to the Bose-Einstein statistics of photons, a classical wave description 
is inadequate to give an exhaustive picture of the phenomenon. In the long term, these gates promise to become 
the basic constituents of quantum computers, either in networked architecture |3, 5], or adopting one-way quantum 
computation approach in order to build cluster states [1, @, H, H, 03 • 

Correct functioning of these devices requires true single photon states. Due to the difficulty of g enerating such 
non classical states [ll|, [l2j , to date parametric down conversion (PDC) has been widely used [H, EE EU , 

where 

high frequency pump photons are non-deterministically converted to pairs of frequency entangled daughter photons 
in a pa ssag e through a non-linear crystal. Independent photons can be produced via two independent PDC processes 
[l7L fl8L [l9]: of the two photons generated in each down conversion event, one is detected to herald the presence of its 
twin which is subsequently sent to the non-classical interferometer. 

It is important to distinguish between imperfections intrinsic to the gate architecture and those which come 
from non-ideality of the soure [l^]. PDC processes is represented generate two photon states in the form 

\ip)— (j + ea\a\ + e 2 (a\) 2 (a\) 2 + o(e 3 )^ |0), with |e| << 1. Higher order terms can not be discriminated without 

photon number resolving detectors, which are not widely available. Crucially, the frequency correlations of the 
interfering photons with their twins can introduce distinguishability and, consequently, welcher weg (which path) 
information. The presence of such knowledge intrinsically reduces the visibility of the interference and the level of 
entanglement of the state generated by the quantu m g ate [20I |2ll [22j • Many efforts have been devoted in designing 
PDC sources with engineered spectral correlations [23. |24|. 



In the present paper, we give a theoretical analysis of the maximal entanglement which can be produced operating 
a linear optical gate with independent PDC photons. This is an important diagnostic for quantum gate engineering, 
in which we wish to know the contribution of the imperfections of the source to the error budget [n|. Note that a 
first detailed sketch of such a calculation can be found in Ref. [IE , where the effect of frequency correlation on 
entanglement swapping and multi-photon entanglement generation is studied. The results presented here represent a 
more general and detailed approach. 



II. NONCLASSICAL INTERFERENCE OF INDEPENDENT PDC PHOTONS. 

Our picture begins with the theoretical treatment of down conversion in Ref. [27| . restricting, for the sake of clarity, to 
the paraxial approximation and considering the pumping with a plane wave beam in mode locking regime. Specifically, 
we want to calculate the optimal visibility in experiments of interference of two single photons, each one coming from 
a frequency entangled pair, given the spectra of the pump fileds. The interferometric apparatus is the one shown 
in Fig. 1: photons 1/ 2, and 3/4 are pairwise produced in two distinct down-conversion processes. Consequently, 
photons 1 and 4 are directly coupled to the detection to provide the trigger signal, while photons 2 and 3 interfere on 
a beam splitter with transmittance T and reflectance R. 
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FIG. 1: (Color online) Schematic Hong-Ou-Mandel interferometer for independent photons. S v {v p ) and S M (/x p ) represent the 
pump spectra for the upper and lower beam respectively. These produce of a photon pair on the modes a\ and a 2 , a\ and a\, 
respectively. Measurement of a photon in the j— th detector is described by the destruction operator Ej. 



Consider, for instance, the PDC event occurring in the upper crystal in Fig.l . The pump beam, traveling on the z 
direction, has a frequency spectrum S v (y p )=s l ,{up)e 1 ^''^ 1 ' 1 '^ around a central value f2 p ; s v {v) describes the distribution 
of the frequencies and 4>y{v p ) their phase relation with respect the central wavelength. In a single event, the pump 
generates of two fields at frequencies lo\ and tu 2 , with wavevectors k\ and k 2 , respectively. In these conditions, the 
two-photon wavefunction is given by, 

\ipu)= J dujpduJi du>2 S(u p - loi - oj 2 )S u (ujp - r2p)a{(wi)a2(cj2)|0), (1) 

up to a normalization factor. Effects due to the finite size of the crystal are neglected. In the other crystal, a 
similar process occur; the second pump beam is chosen to have the same central frequency J7 P , and a spectrum 

The operator describing the photon annihilation due to its measurement is, 



E; 



duj'e 



-4 In 2 



a>Ju% 



(2) 



where j=l,2, 3, 4 and a'j(uj') represents the annihilation operator on mode j at the frequency a/. In the expression 
above we adopt gaussian frequency filters, with central wavelength f2 and full width at half maximum (FWHM) a in 
front of the detectors: these define the coherence time of the photons. We choose the same (degenerate) frequency 
il=fl p /2 and the same FWHM for the four filters, since experiments are usually run near these conditions. Specifically, 
the measured fields at the time t are given by the annihilation operators: 

a[(uj')=e- iu) ' t °ai(uj') 



a 2 (uj')=e- w ' t °(Ta 2 {uj') + Re 1 



a' 3 (uj')=e-^(T*as(uj') 
a' A (u)')=e- iJt °a^'l 



R*e 



'03(0/)) 
- lu ' T a 2 {oj')) 



(3) 
(4) 
(5) 
(6) 



where til=U ~ £i/c, ti optical path of i— th photon, c is the speed of light, and r is the optical delay between the two 
arms [28|. The fourfold coincidence rate is. 



C 4 (t) oc / dt 



J\Ej\i! 12 <g> V>34/ 



(7) 



where dt is shorthand for the differential dt\dt 2 dt^dt A |29j . This integral is taken over the coincidence window we 
choose in our experiment; since typically this is much larger than the coherence time of the photons, and of the 
duration of the pump pulse, we can take integrals on infinite intervals (30j . The calculation is then reduced in 
evaluating the amplitude, 



A(ti,t 2 ,t 3 ,t 4 ;T)-- 



i=i 



34/ ■ 



(8) 
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We use the coordinate transform t\ = t ^' 2 ,t+= t3 ~t, t4, , t v _~t\ — t\,t^_—t\ — t%. We also introduce frequency detuning 
for the PDC: Vi=uJi — f2. Notice that v\ + V2=v p , vj, + v^=^i v . Finally, it is found to be more convenient to express 
the amplitude in terms of the differences V-=V\ — V2, fi-=Vi — V3, thus, 



A(t» + ,t V -,t+,t-;T)= I dv p dfx p d^dfi-S v {iy p )e-(^ypS^ti p )e-(^^'e-(^>- + ^- ) 

\T\ 2 (e-™P t: + e -Vpt% e -iv-t"_/2 e -in-e/2^ _ { e -i^[tl+tl+^^-T) e -i^(tX+tl-^=-+T) x ^ 



e — (*+ — *+H 3 1 ")g- i T-( t ^- t + -l 7- 



We neglect a global phase and normalization constants which are unnecessary for our calculations, since we will deal 
exclusively with ratios. If the integration over the frequencies is carried out we arrive to the expression: 

A{t\X_,t^r_-T) = \T\ 2 (^(^)^(^)e- 5 o(*-) 2 e - 5 o 2 (^) 2 ) + 

- |i?| 2 (F v ((t^ +t1 + - r)/2)F M ((t; + r)/2) e - 5 o(^-*';+ ! ^-) 2 e-^(^-*; + ^+-)^ , 

(10) 

where F v (t) is given by the convolution 

F v {t)= I dt'S^e-^ 1 - 1 ^, (11) 



<5q=<t 2 /32 ln2, and S u (t) is the Fourier transform of the pump spectrum. F p (t) is defined in a similar way for the 
second pump. 

The degree of indistinguishability in HOM interference is measured by its visibility, defined as v—1— C m /CM, with 
C m the coincidence rate at optimal temporal superposition of wavepackets (occurring at a delay To) and Cm the rate 
in absence of interference. We thus can express v in terms of the amplitude A(t v + , t u _, t+, tt; r), 

2\TR\ 2 Re[I s (T )} 
|T| 4 + |i?| 4 I N ' [ ' 

where we have defined, 

I S (r)=§ o J d^dx,^!*^^ (13) 



* N= W J dXl \ Fv ^ 2 j dx2 l F ^ X2 ^ 



(14) 



A physical interpretation of the integral Is( T ) is not straightforward; nevertheless, it is evident that not only the time 
symmetry of each profile is involved, but also non trivial correlations between the two pump pulses. The fact that 
visibility is linked to the pump temporal profile in such a non trivial way, stresses the fact that this interference process 
is much more demanding than the one with correlated photon. The same pulse shape that gives a perfect HOM effect 
with frequency entangled photons, could be only partially suitable to reach high visibilities with independent photons. 
We can derive a simpler form for 2s(0) in the case of infinitely large filters. In this case, we find 



/ s (0)oc J dx 1 \S v {x 1 )\ z \S ll {x 1 )\ z . (15) 
By invoking Parseval's theorem, we can rewrite ()15|) in the frequency domain 

/<j(0)oc [du \Su*S li (u )\*, (16) 



4 




5 10 15 20 25 30 

Filter bandwidth (nm) 



1 




















cr v =3nm 


0.8 












r=3 0.6 
> 0.4 
0.2 














5 


10 


15 


20 


25 30 



Filter bandwidth (nm) 

FIG. 2: (Color online) Visibility vo as a function of filter bandwidth a in absence of phase dispersion. Top: both pump spectra 
are gaussian profiles with the same FWHM a v — cr M . Different plots correspond to different values of <t„. Bottom: the two 
gaussian spectra have different FWHM's. Plots relative to different ratios <j m /ov are shown for the case cr„=3nm. 



where the star denotes the convolution. More generally, a sufficient condition for achieving near perfect visibility is 
that the convolutions have to be solutions of the integral equation: 

F v (x + y)F^x-y)=^ J dz e- u '^ + ^ F u {x + z)F^(x - z). (17) 

It can be demonstrated that the condition of maximum visibility correspond to pump pulses which are much narrower 
in time than the coherence time imposed by the filters. In fact, the solutions of equation (fTTf are F l/ (t)=F^{t)=e~ AS ° t ; 
this are obtained by replacing S(t) with a Dirac distribution. Substitution in (|13D and (1141) leads to the maximum 
visibility v=l and Gaussian modulation of the interference pattern observed with correlated photons [||. The consid- 
erations above indicate which physical process spoils the visibility [25l l2t| . When monochromatic pump beams are 
used, this imposes a very strict correlation to the frequency of the PDC photons. Hence, within the bandwidth of the 
filters some four photons events that do not interfere, since the wavelength of the two trigger photons may be used 
to reconstruct the path information. For instance, consider the event in which, after the beam splitter, photons 1 
and 2 have exactly the frequency f2 p /2, and photons 3 and 4 are slightly detuned at the sides of the filter bandwidth. 
This can be realized by a unique configuration, viz. photons 1 and 2 from the first crystal and photons 3 and 4 from 
the second one. No event with the same wavelengths can be observed when photons 1 and 3 belong to the first pair 
and photons 2 and 4 to the second pair. When the frequency correlation is smoothed by adopting broader pump 
beams, such an information is erased and, consequently high visibility is restored. This descends from the fact that 
we introduced further four photon events which now cancel the ones previously non interfering. This process has been 
individuated as the responsible for degradation of p olarization entanglement in similar experiments of multi-photon 
entanglement creation from entangled pairs [25|, |26| |. 
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As an example, we consider a gaussian spectrum with chirping, 

5 v (^)=e- 41n2 (-) 2 e i (^ ) ^+^ ) ^), (18) 
which gives in the convolution of eq. , 

F u (t)=-===e - +4 *o , (19) 

where ^ v —\ ( 4 j" 2 ~ ^2^) * A similar expression holds for F^(t). Figure 2 show the maximum visibility 

vq= 1 2 y»fl|'] u as a function of the filter bandwidth without chirping. This is plotted for the condition of expected 
maximum visibility tq = 0. We find that high visibility is achieved with narrow filters and short duration hence large 
bandwidth pulses. An imbalance in the bandwidth for the two pulses appears as a further source of distinguishability; 
Fig. 2 also show variations in the visibility curve as the ratio cr^/ay, varies, it is clear that the control on this param- 
eter doesn't need to be extremely accurate to achieve satisfactory values; what primarily influences the visibility is 
the presence of a narrower bandwidth pump pulse, rather than the imbalance of the two pump bandwidths. 

Phase dispersion effects of the second order are far less important. The second order coefficient is responsible 
for a broadening in time of S„(t) with respect the non chirped case. When narrow filters are used, its effects are 
negligible if the ratio (j)^ /(A In 2a~ 2 ) remains of the order of few units. The linear coefficient describes a shift 
on the time axis, and can be use to described time jitter between the two pulses. When the pump beams are 
synchronous {<f>i )= <fix 0) , maximum visibility occurs at r=0; if the following pulses are separated by a jitter time 
Tj (i.e. (j)^=Q,(j)i=Tj) the interferometer is not working at its best delay, consequently, the visibility decreases to 
the value v(Tj). An average over the total measurement time leads to the expected value J dTjv(Tj)P(Tj), being 
P(Tj) a proper distribution for the jitter 



III. POLARIZATION ENTANGLEMENT GENERATION VIA HOM INTERFEROMETRY. 

As said, HOM effect is the basic ingredient for building linear optical gates [3|. In particular, we can adopt 
interference on a partially polarizing beam splitter with \Th\ 2 = 1/3 (|TV| 2 = 1) for the horizontal- H (vertical-F) 
polarization to realize a non-deterministic controlled-sign gate [3, HE Eg] • The model depicted above can be utilized 
to calculate how much frequency entanglement affect the performance of such a gate when it is used for entanglement 
generation. We describe entanglement generation via a partially polarizing beamsplitter with Tjj = 1/3 and Ty = 1, 
where extra losses on vertical modes are induced in order to obtain a correct balancemcnt in the entangled state. In 
this architecture, vertical photons never meet at the surface of the beam splitter, while horizontal photons undergo 
HOM effect. This realizes the controlled interaction, together with post-selection via measurement. If a pair of 
photons with diagonal polarization \D) = (\H) + \V)) impinges on the beam splitter, by post-selecting coincidence 
events we generate the entangled state, 

\<p) = -j=(\VV) + \VH) + \HV)-\HH)), (20) 

since \HH) events are erased by the interference. In the non ideal case, the polarization couples to the wavefunction 
of the pair in the time domain, leading to the following expression for the state of the pair, 

(\A)(\VV) + \HV) + \VH))+e\B)\HH)), (21) 



V3 



c- 



where \A) describes the time behavior without interference, and \B) is the time behavior with interference. Ideally, 
we would have \B) = — \A), and e = 1. The state can be rewritten as, 

-j2= (\A){\VV) + \HV) + \ VH) + ec\HH)) + e^l=W\A^)\HH)) , (22) 

where c = (A\B), and {A 1 -} is an orthogonal vector to \A). By tracing out the temporal part, the polarization density 
matrix is obtained, 

Ptt = ^ + N 2 (l - \c\ 2 ){\HH){HH\) , (23) 



6 



being \ip') the non-normalized vector \VV) + \ HV) + \ VH) + ec\HH). In the case under investigation, we can identify 
up to the normalization, 

<0|A) =^(i;)F M (^)e- 5 o( t -) 2 e - 5 o(^) 2 ; (24) 
(0\B) = A(t» + X-,t»+X-;T). (25) 

Consequently, we find, 

e = ^5(1 - v), (27) 

in the case in which Is (r) is a real function. It is straightforward to check that c and e reach their ideal values in cor- 
respondence of the optimal visibility ^=0.8. A simple relation between the visibility and the degree of entanglement 
in the state (j2"3"j) . quantified via the tangle [3l| is found, 
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FIG. 3: (Color online) Plot of the tangle in the polarization state of eq. (|23[) as a function of the HOM visibility, normalized 
to its maximum value vm = 0.8 



The plot in Fig. 3 shows the behavior of this curve as a function of the relative visibility v/vid; it is clear that the 
degree of entanglement rapidly is degraded as the visibility decreases. Tangle remains well under the value 0.9, even 
for values of the visibility which are commonly considered more than satisfactory in experiments. Thus, the maximum 
attainable entanglement strongly depends on the choice of the pump and the detection bandwidths, as shown in Fig. 
4 for the case of Gaussian pumps. A comparison of the figures indicates that tangle is a much stricter criterion in 
designing an experiment. 



IV. DEPARTURES FROM THE IDEAL CASE. 



The results obtained above neglect the finite size of the crystal and the effect of phase dispersion of the PDC 
photons; therefore they remain valid in the limit of short crystal. 

So far we have disregarded the fact that the crystals have finite thickness. If this is taken into account, the state 
([1]) has to be modified as 

1^12)= / dujpduji dui 2 S(ujp —ui — uj 2 )A„(ujp - Q p ,oji - uj 2 )S 1/ (uj p - fl p )a\ (u>i)al(u) 2 )\0) , (29) 
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FIG. 4: (Color online) Comparison of the plots of the relative visibility and the tangle in the polarization state in eq. (|23[) 
as a function of the pump bandwidth in the case of a Gaussian pump. Different curves correspond to different values of the 
detection bandwidth. 



where A. u (i/ p ,i/^) is the function describing non-perfect phase matching. The amplitude (|10p is modified by substi- 
tuting F^(t)e _4t5 o* 2 with the convolution, 

G„(t,t / ) = I ' dsds'F v (s)e- i5 o s " 'A„(s-t,s' -if), (30) 



where A„(t, t') is the two dimensional Fourier transform of A„(i/ p ,i/_). A similar substitution should be carried out 
for the second pump beam. 

In the case of degenerate type II phase matching, this can be approximated as [27j |. 

A„(i/p, v-) = sinc(a p ^ p + a_i/_), (31) 

being a p and a_ constants which depend on the length L and the dispersion in the crystal, viz. 

1 1 



4 \u (n) u e {n)J 2 Ue {n p y 

L ( 1 1 ( "~' 



Oi- = — 



2 \u (ty U e (fi) J ' 

where u a (Q) ( u e (f2)) is the group velocity dispersion for an ordinary (extraordinary) beam at the frequency Cl. In 
this case, the expression for its transform is, 

1 t a 

A„(t, t') = — rect(— )S(t' - — t), (33) 
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which leads to the result, 

G„(i,f) = — I"* 12 dsF L ,(t + s)e~ 4S "( t ' + ^ s y . (34) 

a P J-a p /2 

Note that in the short crystal limit a p — > 0, we restore the expression G v (t, t') — I*i,(i)e — 4 *°(* ' 2 . In the opposite limit 
of infinite crystal the correlation (|3Tj) becomes perfect 

v v = (35) 

a p 

Being a p , a_ linearly dependent on the crystal length L, the ratio — remains finite when passing to this limit. Such 
a correlation implies that a PDC pair with a frequency mismatch v_ can be generated only by the component of the 
pump beam detuned of ——v p . Therefore, it can reduce the advantage of using a broad spectrum pump, since the 
crystal itself limit the matching bandwidth. As an example, in Fig. 5 the results for degenerate type II PDC from 
a 400nm pump pulse are shown for two widely adopted materials: LiB 3 Os (LBO), and BaB 2 04 (BBO). The pump 
bandwidth is 5nm, while the detection filters are chosen to have 2nm FWHM: in the ideal case, these values give 
near perfect visibility. It is clearly shown that the finite size of the crystal can be neglected for thicknesses in the 
range 2mm to 4mm. This is the size of bulk crystal slabs commonly adopted in experiments. At these wavelength, 
the adoption of LBO instead of BBO reduces the distinguishability effect. Anyway that would cause a significant 
reduction in the brightness of the source, given the smaller nonlinearity of LBO. 

As a final remark, we note that when using type II PDC it is convenient to use the signal photons as triggers and 
the idler photons in the interferometer. In fact, distinguishability of signal and idler photons can arise from different 
dispersion in the PDC crystal due to birefringence. 
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FIG. 5: (Color online) Plot of maximum achievable visibility as a function of crystal length in Type II degenerate 400nm^800nm 
PDC, for two widely adopted materials, LBO and BBO. The pump bandwidth is 5nm, the detection bandwidth is 2nm. This 
choice gives near unity visibility in the ideal case. Sellmeier coefficients are taken from Ref. [32] 



V. CONCLUSION 



The results we have presented indicate that when running entangling gates with independent PDC photons, large 
bandwidth pulses are somehow preferable to narrow ones. Their adoption reduces the welcher weg knowledge in the 
interference process coming from frequency entanglement. We showed that visibility has to be considered a forgetful 
parameter in the design of these experiments. The non trivial relation between tangle and visibility implies that the 
conditions to fulfill to attain a high level of entanglement are stricter than those for reaching near perfect visibilities. 
Such an effect of frequency correlation can be easily reduced by properly choosing pumping and detection conditions 
with existing materials. Further investigation of this process can include the effects of group velocity dispersion on 
the arrival time of the PDC photons, resulting in a time jitter [HJ 

At a first sight, it might seem counterintuitive that broadband pumping is better than single mode pumping; in 
fact it has been commonly observed that interference of dependent photons in pulsed regime is more difficult to reach 
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than with cw pumping. Nevertheless, using a mode locked pump per se does not present in principle limitations. 
For the independent photon case, the actual temporal shape of the pump pulse is crucial, since knowledge comes 
from non trivial correlations between the pump pulses. In real life implementations, other sources of imperfections 
may impose restrictions on pump and detection bandwidth, e.g. dispersion in dielectric mirrors, wavelength-sensitive 
performance of beam splitters. Furthermore, restricting the duration of the pump pulses is less resilient to jitter and 
narrow detection filter reduce the apparent brightness of the source, leading to longer measurement times. 
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Note added 

During the refereeing process of this manuscript, we become aware of a paper presenting a similar treatment (34j . 
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